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SUMMARY 


This report documents the derivation and definition of a linear aircraft model for a rigid aircraft of constant 
mass flying over a flat, nonrotating earth. The derivation makes no assumptions of reference trajectory or 
vehicle symmetry. The linear system equations arc derived and evaluated along a general trajectory and 
include both aircraft dynamics and observation variables. 

INTRODUCTION 

The need for linear models of aircraft for the analysis of vehicle dynamics and control law design is well 
known. These models are widely used, not only for computer applications but also for quick approximations 
and desk calculations. Whereas the use of these models is well understood and well documented, their 
derivation is not. The lack of documentation and, occasionally, understanding of the derivation of linear 
models is a hindrance to communication, training, and application. 

This report details the development of the linear model of a rigid aircraft of constant mass, flying over a 
flat, nonrotating earth. This model consists of a state equation and an observation (or measurement) equa- 
tion. The system equations have been broadly formulated to accommodate a wide variety of applications. 
The linear state equation is derived from the nonlinear six-degrcc-of-freedom equations of motion. The 
linear observation equation is derived from a collection of nonlinear equations representing state variables, 
time derivatives of state variables, control inputs, and flightpath, air data, and other parameters. The linear 
model is developed about a nominal trajectory that is general. 

Whereas it is common to assume symmetric aerodynamics and mass distribution, or a straight and level 
trajectory, or both (Clancy, 1975; Dommasch and others, 1967; Etkin, 1972; McRuer and others, 1973; 
Northrop Aircraft, 1952; Tlielander, 1965), these assumptions limit the generality of the linear model. The 
principal contribution of this report is a solution of the general problem of deriving a linear model of a rigid 
aircraft without making these simplifying assumptions. By defining the initial conditions (of the nominal 
trajectory) for straight and level flight and setting the asymmetric aerodynamic and inertia terms to zero, 
one can easily obtain the more traditional linear models from the linear model derived in this report. 

Another significant contribution of this report is the derivation and definition of a linear observation 
(measurement) model. The observation model is often entirely neglected in standard texts. A thorough 
treatment of common aircraft measurements is presented by Gainer and IIolTman (1972), and Graccy (1980) 
provides a detailed discussion of speed and altitude measurements. However, neither of these references 
present linear models of these measurements. This report relics heavily on these two references and uses their 
results as one of the bases for the nonlinear measurement equations from which the linear measurement 
model is derived. Also included in this report is a large number of other measurements or variables for 
observation that have been found to be useful in vehicle analysis and control law design. 

Duke and others (1987) describe a FORTRAN program called LINEAR that derives a linear aircraft 
model by numerical differencing (Dieudonne, 1978). The program LINEAR produces a linear aircraft model 
(both state and observation matrices) that is equivalent to the linear models defined in this report. 

This report is divided into two main sections that define the reference systems and nonlinear state and 
observation equations (section 1) and derive a linear model presented in the appendixes (section 2). The 
appendixes contain a definition of the linear aerodynamic model used in this report (app. A), a derivation 
of the wind axis translational acceleration parameters (app. B), generalized linear derivatives of the non- 
linear state and observation equations (app. C), and the individual derivatives of the stale and observation 
equations (app. D). The details of the principal results of this report arc presented in appendix D. 



SYMBOLS 


A total aerodynamic axial force, lb 

a speed of sound, ft/sec 

« n normal accelerometer output, g 

a u ,i output of normal accelerometer not at vehicle center of gravity, g 

a x output of accelerometer aligned with vehicle body x axis, g 

ci Xt i output of accelerometer aligned with body x axis, not at vehicle center of gravity, g 

a Xi k kinematic acceleration in vehicle body x axis, g 

c i y output of accelerometer aligned with vehicle body y axis, g 

a y> i output of accelerometer aligned with body y axis, not at vehicle center of gravity, g 

ciy^ kinematic acceleration in the vehicle body y axis, g 

a z output of accelerometer aligned with vehicle body z axis, g 

(i z ,i output of accelerometer aligned with body z axis, not at vehicle center of gravity, g 

kinematic acceleration in vehicle body z axis, g 
b reference span, ft 

C \ generalized force or moment coefficient 

C^ x derivative of generalized force or moment coefficient with respect to arbitrary variable x 

c reference aerodynamic chord, ft 

D total aerodynamic drag, lb 

D X I Z -Iy 

Dy I x ~ I z 

D Z Iy ~ I X 

E s specific energy, ft 

F arbitrary force or moment 

fpa flight path acceleration, g 

g acceleration due to gravity, ft /sec 2 

go acceleration due to gravity at sea level, ft/sec 2 

h altitude, ft 

h,» altitude measurement not at vehicle center of gravity, ft 

I inertia tensor 

I x moment of inertia about x body axis, slug-ft 2 

I xy product of inertia in x-y body axis plane, slug-ft 2 

Ixz product of inertia in x-z body axis plane, slug-ft 2 

I y moment of inertia about y body axis, slug-ft 2 

lyz product of inertia in y-z body axis plane, slug-ft 2 

I z moment of inertia about z body axis, slug-ft 2 

h 

^2 Ixy^z “I" lyz Ixz 

I3 Ixy lyz H - lylxz 

h hh - 4 

lb Ixlyz 4 " Ixy Ixz 

h U y ~ Ily 

L total moment about 

t unit length, ft 

M total moment about 


a: body axis, ft-lb; or, total aerodynamic lift, lb 
y body axis, ft-lb; or, Mach number 


2 


m 

N 

n 

Ps 

p 

Pa 

Ps 

Pi 

7 

7 

7c 

7c/ Pa 
7s 
Re 
Re' 


r s 
S 
T 
T t 

t 

u 

V 

v 

w 

X a 

X T 

X 

Y 

la 

Xg 

It 

y 

Z a 

Zt 

z 

a 

a,i 

P 

P,i 

1 

Si 

e 

/* 

p 

<& 


vehicle mass, slugs 

total moment about 2 body axis, ft-lb; or, total aerodynamic normal force, lb 

load factor 

specific power, ft/sec 

roll rate (about x body axis), rad/sec 

static or free-stream pressure, lb/ft 2 

stability axis roll rate, rad/sec 

total pressure, lb/ft 2 

pitch rate (about y body axis), rad/sec 

dynamic pressure, lb/ft 2 

impact pressure, lb/ft 2 

Mach meter calibration ratio 

stability axis pitch rate, rad/sec 

Reynolds number 

Reynolds number per unit length, ft 
yaw rate (about z body axis), rad/sec 
stability axis yaw rate, rad/sec 
surface area of wing, ft 2 

total angular momentum; or, ambient or free-stream temperature, R 

total temperature, °R 

time 

velocity along x body axis, ft/sec 
vehicle velocity, ft/sec 
velocity along y body axis, ft/sec 
velocity along z body axis, ft/sec 
total aerodynamic force along x body axis, lb 
total gravitational force along x body axis, lb 
total thrust force along x body axis, lb 
vehicle position along x earth axis, ft 
total aerodynamic sideforce, lb 
total aerodynamic force along y body axis, lb 
total gravitational force along y body axis, lb 
total thrust force along y body axis, lb 
vehicle position along y earth axis, ft 
total aerodynamic force along z body axis, lb 
total gravitational force along z body axis, lb 
total thrust force along z body axis, lb 
vehicle position along z earth axis, ft 
angle of attack, rad 

angle-of-attack measurement not at vehicle center of gravity, rad 
angle of sideslip, rad 

angle-of-sideslip measurement not at vehicle center of gravity, rad 

flightpath angle, rad 

ith control surface deflection 

pitch angle, rad 

coefficient of viscosity, lb/ft-sec 
density of air, lb/ft 3 
arbitrary function 



<j> bank angle, rad 

rj) heading angle, rad 

Vectors 

a 
E 
F 
f 

g 
H 
h 
M 
R 
u 
V 
x 

y 

8 u 
<5x 
Sic 

Matrices 

A state matrix of the generalized state equation, Cx = /lx + 13 u 

A f state matrix of the state equation, x = A'x -f B'u 

B control matrix of the generalized state equation, Cx = Ax + B u 

B' control matrix of the state equation, x = A'x + B'u 

C system matrix of the generalized state equation, Cx = Ax + B u 

F feedforward matrix of the generalized observation equation, y = II x + Cx + Fu 

F' feedforward matrix of the observation equation, y — IFx + F'u 

G derivative observation matrix of the generalized observation equation, y = II x -f Gx + Fu 

II observation matrix of the generalized observation equation, y — II x + Gx + Fu 

IF observation matrix of the observation equation, y = II'x + F'u 

/ intertia tensor 

J ' scaling matrix for inertia tensor 

Fbv transformation matrix from earth to body axes 

R transformation matrix from earth to body axes 

T angular velocity matrix in the generalized state equation, Tx = f[x(Z), x(2), u(<)] 

On x ? 7 i 71 x m matrix of 0 values 

1 7 i x 7 n an n X m matrix with values of 1 on the diagonal 


body axis acceleration vector 

attitude vector of Euler angles 

total force vector 

state vector function 

observation vector function 

total angular momentum vector 

sum of higher order terms in Taylor series 

total moment vector 

position vector in earth axis system 

input or control vector 

vehicle velocity vector 

state vector 

observation vector 

perturbation of control vector 

perturbation of state vector 

perturbation of time derivative of state vector 

rotational velocity vector 


Subscripts 


a 

b 

D 

g 

h 

h 

,i 

,k 

L 

l 

m 

n 

n 

P 

s 

T 

t 

V 

w 

X 

xy 

xz 

Y 
y 

yz 

z 

0 


aerodynamic; or static or, free stream 

body axis system 

drag 

gravitational 

displacement of altitude instrument 
displacement of altitude rate instrument 
not at vehicle center of gravity 
kinematic 
lift 

rolling moment 

pitching moment 

yawing moment 

orthogonal 

power plant induced 

stability axis; or, specific 

thrust 

total 

vehicle-carried vertical axis system 
wind reference axis system 
displacement in x body axis 
x-y body axis plane 
x-z body axis plane 
sidcforce 

displacement in y body axis 
y-z body axis plane 
displacement in the z body axis 

at sea level, standard day conditions; or, nominal conditions 


Superscript 

T transpose 

1 NONLINEAR SYSTEM EQUATIONS 

The motion of an aircraft as a rigid body can be described by a set of six nonlinear simultaneous second- 
order differential equations. These equations, representing the translational and rotational motion of the 
vehicle, can be formulated in the notation of Kwakernaak and Si van (1972) and Dieudonne (1978) as a 
time-invariant system expressed as 

x(f) = f[x(f), u(<)] (1-1) 

where x(Z) is the 12-dimensional time- varying state vector ( t being time), x(/) is the derivative of x(/) with 
^respect to time, u (t) is the /c-dimensional time-varying input or control vector, and f is a 12-dimensional 
nonlinear function expressing the six-degree-of-freedom rigid body equations. 

Measurements of the vehicle state can be represented by the observation equation 

( 1 - 2 ) 
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y(0 = g[x(/),u(0] 



where y(t) is an ^-dimensional time- varying observation vector and g is an £- dimensional nonlinear func- 
tion expressing the relationship of the true vehicle state and control vectors to the observed parameters. 
Typically, the function g characterizes the dynamics and location of the sensors. 

For the aircraft analysis and design problem, both the nonlinear and linear system equations are formu- 
lated more broadly than just described (Edwards, 1976; Maine and Iliff, 1980, 1986). The nonlinear system 
equations include x(/) terms in both the state and observation functions. In fact, in the most extended 
form the state equation is expressed in terms of transformed variables (discussed in section 1.2.1). These 
generalized equations form the basis of the analysis in this report. The generalized system equations arc 


Tk(t) = f[x(<),x(<),u(<)] 

(1-3) 

y (0 = gM*),x( 0 ,u(z)] 

(1-0 


where T is a constant 12 X 12 angular velocity matrix. 

1.1 Definition of Reference Systems 

While numerous reference systems are used in aerospace applications, this report is limited to four reference 
systems: the body, the wind, the vehicle-carried vertical, and the topodetic reference systems. The stability 
axes are also defined even though this reference system is used only to define the stability axis rotational 
rates (section 1.3.8). 

Within this report the translational equations are referenced to the wind axes, and the rotational 
equations are referenced to the body axes. Measurement equations are primarily referenced to the body 
axes when the use of a reference system is needed. The use of this mixed axis system definition in both 
the nonlinear and linear models is related to the measurability and meaningfulness of quantities. Because 
the aerodynamic forces act in the wind axes, this reference system is used for the translational equations. 
For instance, angle of attack, velocity, and angle of sideslip are either directly measurable or closely related 
to directly measurable quantities, while the body axis velocities (u, u, and w in the and 2 directions, 
respectively) are not. The body axis rotational rates are measured by sensors fixed in the body axes; wind 
axis rates can be derived only from these quantities through axis transformations. 

The first reference system to be described is the topodetic reference system, also called the earth-fixed 
reference frame (Etkin, 1972), the earth axes (Thelander, 1965), and the Eulerian axes (Northrop Aircraft, 
1952). The topodetic reference frame is considered fixed in space (and hence, inertial) with the orientation 
of the axes as shown in figure 1; the x axis is directed north, the y axis east, and the z axis down. The 
vehicle position (x and y) and altitude (/i) are measured from the origin of this reference system. 

The vehicle-carried vertical axis system (fig. 2; Etkin, 1972) has its origin at the center of gravity of the 
vehicle. The x v axis is directed north, the y v axis east, and the z v axis down. This axis system is obtained 
by a translation of the topodetic axis system to the vehicle center of gravity. The attitude of the aircraft 
(heading, pitch, and bank angles 0, and <^>, respectively) is described in terms of the orientation of the 
aircraft body axes with respect to the vehicle-carried vertical axes. 

The origin of the body axis system (fig. 3) is the vehicle center of gravity. The x axis is directed toward 
the nose of the aircraft, the y axis toward the right wing, and the z axis toward the bottom of the aircraft. 
The specific orientation of the actual body axes relative to the vehicle body is somewhat arbitrary. For 
symmetrical aircraft, the .t and ^ axes are in the plane of symmetry; for asymmetrical aircraft, these axes 
are located in a plane approximating what would be the plane of symmetry. The positive direction for the 
body axis rates (roll, pitch, and yaw rates, p, q y and r, respectively), the body axis velocities (u, and w ), 
and the body axis moments (Z, A/, and N about the ?/, and z axes, respectively) are shown in figure 3. 
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Figure 1. Topodetic axis system. Figure 2. Relationship between topode- 

tic and vehicle-carried vertical axis sys- 
tems , 



♦ 
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Figure 3. Body axis system . 
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The relationship between the vehicle- carried vertical and body axes is shown in figure 4. The Euler 
angles (^, 8, and <p) define the orientation of the body axes with respect to the vehicle-carried vertical 
axes. The rotations required to transform the vehicle-carried vertical axes to the body axes are shown in 
figure 5. The heading angle ip is a rotation about the z vehicle- carried vertical axis into a new axis system 
(designated (.ti, ?/i, ^x) in fig. 5); the pitch attitude 6 is a rotation about the yi axis into the (# 2 , 1 J 2 , ^ 2 ) 
axes system; the roll attitude <p is a rotation about the 1 J 2 axis into the body axes. 


x 2 , x b 



Figure 4- Relationship between vehicle- carried vertical and body axis systems. 


These rotations arc described by 






L<p 


and the total rotation is described by 


COS Ip - 

- sin ip 0 

sin ip 

cos ip 0 

0 

0 1 

cos 8 

0 sin 8 ' 

0 

1 0 

— sin 9 

0 cos 8 

1 0 

0 

0 cos <p 

— sin <P 

0 sin (p 

cos <p 


— L yj L@ L £ — 


cos 8 cos ip 
sin <p sin 8 cos ip 
— cos (p sin ip 
cos sin 9 cos ip 
+ sin 4> sin ip 


cos 9 sin ip — sin 8 
sin <p sin 8 sin ip sin <p cos 8 
+ cos <p cos ip 

cos <p sin 8 sin ip sin <p cos 8 
- sin (pc os ip 


(1-5) 


( 1 - 6 ) 


(1-7) 


( 1 - 8 ) 


s 



Because Xgv is a unitary matrix, the transformation from the body axes to the vehicle-carried vertical 
axes is ^BV‘ 

The relationships between the body, wind, and stability axes are shown in figure 6. All three axis 
systems have their origin at the center of gravity of the aircraft. The x axis in the wind reference system 
(x w ) is aligned with the velocity vector of the aircraft. The angle of sideslip (3 and angle of attack a define 
the orientation of the wind axes with respect to the body axes. (The stability axes are shown in figure 6 
also. This reference system is displaced from the wind axis system by a rotation / 3 and from the body axis 
system by a rotation —a.) 



Figure 6 . Relationship of body, stability, and wind axes . 

Also shown in figure 6 are the components of the velocity vector V in the body axes ( u , v , and w) 
and the definition of positive rotations for a and /?. It should be noted that j3 is a positive rotation in a 
left-handed coordinate system, whereas the positive sense of all other rotations used in aircraft analysis are 
positive in a right-handed coordinate system. 

The definitions of the body axis velocities (fig. 6) are 


The total velocity V, angle of attack a, and 
velocities as 


u = V cos a cos (3 

(1-9) 

v = V sin (3 

(MO) 

w = V sin a cos /? 

(1-11) 


angle of sideslip (3 can be expressed in terms of these body axis 


V 


a 


P 


|V| = 0 2 + V 2 + w 2 ) 1 / 2 


-l 


tan 


sin 1 — 


w 

u 

v 

V 


( 1 - 12 ) 

(1-13) 

(1-14) 
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1.2 Nonlinear State Equations 


For the aircraft problem, the state vector x is 12 x 1 vector composed of four 3x1 subvectors representing the 
vehicle rotational velocity, the vehicle translational velocity, the vehicle attitude, and the vehicle location: 


= [xM x? xJP 

(1-15) 

Xi = [p q r] T 

(1-16) 

x 2 = [V a /?] T 

(1-17) 

x 3 = [<t> 0 0] T 

(1-18) 

X 4 = [h x y] r 

(1-19) 


with xi, x- 2 , x 3 , and x 4 being the rotational velocity, translational velocity, attitude, and position subvectors, 
respectively. The vehicle rotational and translational velocity are defined within the aircraft-fixed axis 
systems. In the formulation of the state used in this report, the vehicle rotations are body axis rates, whereas 
the vehicle velocity terms are stability axis parameters. The vehicle attitude and location parameters are 
earth relative. 

The vector function f, relating the state vector its time derivative, and the control vector to the time 
derivative of the state vector with respect to time, is a 12-dimensional vector function composed of four 
3-dimensional vector subfunctions: 


f[x(t),k(t)Mt)} = [f7fJfJfI] T d-20) 

where f 1? f 2 , f 3 , and f 4 are the vector functions that relate the x(/), x(/), and u(/) vectors to the rotational 
acceleration, translational acceleration, attitude rate, and earth-relative velocity subvectors of x(/). In the 
following sections, each of these subfunctions will be developed separately. The details of the derivation 
of these subfunctions can be found in any of the standard references on aircraft dynamics (Etkin, 1972; 
Meitner and others, 1973; Thelander, 1965). 


1.2.1 Rotational acceleration. — The subfunction fi of f from which the rotational acceleration 
terms in the x vector are derived is based on the moment equation 


M=4-H (1-21) 

at 

where M is the total moment on the vehicle and H is the total angular momentum of the vehicle. This 
expression can be expanded to 

M=f(/fi) + fix(/fi) (1-22) 

St 

where 8 /6t is the time derivative operator in a moving reference frame (such as the vehicle body axis system) 
and the substitution 

H = /ft (1-23) 

has been used to replace the total angular momentum term with the product of the inertia tensor I and 
the rotational velocity vector ft. (The inertia tensor is assumed to be constant with time.) The definition 
of the terms in equation (1-22) follow: 





L + L t ' 

M = 

ZM 

= 

M + Mj 


ZN 


N + Nj 


(1-24) 
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with Z, M, and iV being the aerodynamic total moments about the x, ?/, and 2 body axes, respectively, and 
L t, and A r x the sums of all power-plant-induced moments; 

lx Ixy 

I — ~Ix y ly 

— Ixz lyz 

where J x , 7 y , and I z are the moments of inertia about the x, y, and z body axes, respectively, and I xy , l xz% 
and I yz are the products of inertia in the x-?/, x-z, and y-z body axis planes, respectively; and 

fl = Xi = [p q r] T (1-26) 

where p, q , and r are the rotational rates about the x, y , and z body axes, respectively. Because it is 
assumed that the inertia tensor is a constant with respect to time, equation (1-22) can be rewritten as 

7 fi = /- 1 ( M-flx/fi) (1-27) 


-h 

-I, 


yz 


(1-25) 


This is the vector subfunction for the rotational acceleration. Designating this subfunction as f l7 the 
following definition applies: 

f i[x(0,x(0iU(0] = r _ 1 [M- ft X (Til)] (1-28) 

where 


j- t Q = { i[ X (<),x(0,u(0] 


(1-29) 

(1-30) 


Since the inverse of the inertia tensor / Ms given by 


I- 1 


1 

det I 


h h h 

1 7 1 4 hi 

h h h 


(1-31) 


where 


det I = 

h = 
h = 
h = 
u = 
h = 
h = 


I x Ty I z Ixlyz Izlxy ~ ly^xz ~ 2 Tyzlxzlxy 

(1-32) 

I y I Z ~ Iyz 

(1-33) 

Ixylz “1" lyzlxz 

(1-34) 

Ixy Iyz “1" lylxz 

(1-35) 

1 I - I 2 

l x 1 z l xz 

(1-36) 

Ixlyz ~\~ Ixylxz 

(1-37) 

I J _ /2 

1 !/ 

(1-38) 
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the expression for the rotational accelerations can be expanded as a set of scalar equations. 

P = rt-p Lh + SM/i + S Nh - P 2 {lxzh - Ixyh ) + P'lihzh ~ Iyzh ~ D Z I 3 ) 
det I 

— pr(I xy Ii + Dyh — lyzl 3 ) + 1 ~ 4y4) ~ < l r {D x I\ — Ixyh + hzh) 

— r 2 (I yz Ii - I XZ I 2 )] 0'39) 

q = [E LI 2 + SM/ 4 + EA r J 5 - p 2 (4*4 - Jxyis) + P9(^*/2 - Iyzh ~ D *h) 

det I 

— pr(I X yh + D y L 1 — Iyzh ) + q 2 (Iyzh ~ Ixyh) — qr(D x I 2 — Ixyh + hzh ) 

— r\lyzh - 4*4)] (!- 4 °) 

f = — — — -[SZ /3 + SM/s + SiV / 6 - p 2 (4*4 - Ixyh) + P<7(4*4 - lyz h ~ D z I g ) 
det / 

— P r (Ixyh + Dyh ~ Iyzh) + q 2 {Iyzh ~ Ixyh) ~ Q r (D x h — Ixyh + Ixzh) 
-T 2 (I yz h-Ixzh)} C 1 ' 41 ) 

where 

D x = I z -Iy (1-12) 

D y = I X ~ I Z (1-13) 

D z = Iy- I x (1-44) 


Equation (1-3) defines the generalized nonlinear state equations as 

Tk{t) - f[x(/),x(<),u(0] 

This equation, although more complicated than the nonlinear equations defined by equation (1-1), allows for 
a more tractable formulation of the state equation by using the matrix T to provide a means of addressing 
the rotational accelerations in a decoupled axis system. 

The derivation of the rotational acceleration terms is based on the moment equation (1-22): 


C 

m = —(/ft) + ft x m 

ot 

Rearranging terms and assuming that the inertia tensor is constant with respect to time, the equation can 
be written as 

/|-ft = M-ftx/ft (1-45) 

ot 

The rows of this vector equation are now scaled using the following scaling matrix: 


J' = 


1/4 0 0 
o i// y o 
o o 1/4 


(1-46) 


This matrix, when premultiplying equation (1-27), merely divides the first row by the roll inertia I x , the 
second row by the pitch inertia I y , and the third row by the yaw inertia I z . Using tire definition 

J = J'l (1-47) 
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the resulting equation is 


(1-48) 


J-Q = J'M - J'(Q x IQ) 
ct 


and J can be written as 


J = 


1.0 

Ixy! Iy 
hz/Iz 


^xyj^x 

1.0 

lyzj 


Jxzl lx 

~ lyzl ly 

1.0 


Equation (1-48) can be expanded and expressed as 


p' 


i 

i 

o 

Lfi 


P 


=r 

Ixy/Iy 1-0 Iy Z /Iy 


<i 

r' 


. -Ixzlh -Iyz/Iz 1.0 


f 


XL/lx r plxy/lx 4" PQIxz/I : r H" TQly/ lx l ^)lyz/^x Q^Iz/^x 

TiM/Iy - rplxfly + rql xy /ly - pql yz /ly + (r 2 - p 2 )I xz /I y + prl z /l y 
E Nfl z + qplx/lz - qrlxz/Iz + prl yz /l z + (p 2 - q 2 )I xy jl z - pqly/h 

where p\ </', and r f are the decoupled rotational accelerations of the vehicle. 

Using the definition of J in equation (1-49), the matrix transformation T can be defined as 


(1-49) 


(1-50) 



J 1 03x3 
1 

1 

1 

1 

06X6 

! 

03X3 | 13x3 

1 

1 

_ I 

to 

X 

CO 

O 

1 

1 

1 

1 

1 

1 

lcx6 


(1-51) 


which would be an identity matrix except for the presence of the inertia terms in the upper left-hand corner. 
Thus, the vector subfunctions for the generalized state equation defining vehicle translational acceleration, 
vehicle attitude rates, and earth- relative velocities are the same as those defined for the standard nonlinear 
state equations in sections 1.2.2, 1.2.3, and 1.2.4, respectively. 


1.2.2 Translational acceleration. — Derivation of the translational acceleration vector subfunction 
f 2 is based on the force equation 

F= ^ (mV) (1*52) 

where F is the total force acting on the vehicle and m is the vehicle mass. This expression can be expanded to 

F = m (*V + ft X v) (1-53) 

with the assumption of constant mass with respect to time and the following definitions of F and V: 

F = [EX E Y EZ] T (1-54) 

where ELY, E Y, and E Z are the sums of the aerodynamic, thrust, and gravitational forces in the x , y, and 

2 body axes, respectively, and 

V = [tx V u>] T (1-55) 
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(1-56) 


Rearranging the terms of equation (1-52) gives an expression for the translational acceleration: 

l v = -f - n x v 

St nr 

This equation expresses body axis accelerations in terms of body axis forces, angular rates, and velocities. 
However, the desired form of this relation requires the translational accelerations in the wind axis system; 
that is, in terms of the magnitude of the total vehicle velocity V , angle of attack a, and angle of sideslip /?, 
which are expressed by equations (T9) to (Til) 

u = V cos a cos /3 
v = V sin P 
w — V sin q cos (3 


and equations (1-12) to (1-14) 


V = \V\- ( u 2 + V 2 + w 2 ) 1/2 



P 



The wind axis translational acceleration terms (derived in app. B) are summarized as: 

[l>d/j] T = f 2 [x(0,*(0,u(0] (1-57) 

where 

V — — [ — D cos P + Y sin /? + A^x cos a cos p + Yj sin P + Zx sin a cos /? 
m 

— mg(cos a cos P sin 9 - sin P sin <p cos 6 - sin a cos P cos $ cos 0)] (1-58) 

a = -[ — L + Zx cos a — Ax sin a -f mucosa cos<^cos# + sin a sin 9)] 

Vm cos P 

+ q - tan /3 (pcos a + r sin a) (T59) 

p = —^—[1? sin P + y cos p — At cos q sin P + lx cos P - Zx sin a sin P 
mV 

+ mg(cos a sin P sin 0 + cos P sin <f> cos 9 — sin a sin p cos ^ cos 0)] -f p sin a — r cos a (1-60) 

with D being total aerodynamic drag; Y total aerodynamic sideforce; and Ax, Tt? an( J Zx total thrust 

force along the x , ?/, and z body axes, respectively. 


1,2.3 Attitude rates. — The matrix R that transforms angular velocities in the earth-fixed axis system 
into body axis angular velocities is defined by 


R = 


1 0 -sin# 

0 cos <i> sin $ cos 9 

0 — sin (f> cos 4> cos 9 


( 1 - 6 - 1 ) 


lb 



where R is derived by Maine and TlifF (1986) from the total angular velocity of the aircraft expressed in 
terms of the derivatives with respect to time of the Euler angles (<£, 0, ip): 


V 


4> 


'i 

0 0 


'o' 


'l 

0 0 


cos 9 0 

— sin 8 


'o' 

9 

— 

0 

+ 

0 

cos (p sin (p 


6 

+ 

0 

cos (p sin <p 


0 1 

0 


0 

r 


0 


0 

— sin (p cos (p 


0 


0 

— sin <p cos (p 


sin 9 0 

cos 9 




1 0 — sin 9 


4> 

0 cos <p sin <p cos 9 


6 

0 — sin cp cos (p cos 9 




This transformation from earth-fixed to body axes can be expressed by the equation 


(1-62) 


q - r (i E ) 


(1-63) 


where E is an attitude vector whose components are the Euler angles: 


e = [<f) e xi>] T 


(1-64) 


Premultiplying both sides of equation (1-63) by R 1 and rearranging terms yields the equation for the 
attitude rates, 

(1-65) 

which can be expanded into the scalar equations 

4> — V + <7 sin <p tan 9 + r cos <p tan 5 (1-66) 

6 = q cos 4> — r sin <p ( 1-67) 

jp = qsin(p sec 9 + r cos <p sec 6 (1-68) 


1.2.4 Earth-relative velocity. — The matrix T B v that transforms earth axis system vectors into the 


T'BV 


cos ip — sin ip 0 
sin ip cos ip 0 

0 0 1 


as 



cos 9 

0 

sin 9 

0 

1 

0 

— sin 9 

0 

cos 9 


0 

— sin (p 
cos <p 


1 0 
0 cos <p 
0 sin <p 

cos 9 cos ip cos 9 sin \p — sin 9 

sin <p sin 9 cos ip — cos <p sin ip sin <p sin 9 sin ip + cos <p cos ip sin <p cos 9 
cos (p sin 9 cos ip + sin (p sin ip cos <p sin 9 sin ip - sin <p cos ip cos<p cos 9 

The specific relationship between earth- relative velocities and body axis velocities is expressed by 


(- 

\dt 


r ) 


V - Z B v 

where R is the earth axis system vector defining the location of the vehicle: 

R = h z] T 


(1-69) 
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(1-70) 



(1-71) 


with z — —h. 

The equation for the earth-relative velocity can be formulated as 





in which these velocities are expressed in terms of body axis velocities. Using equation (1-72) and the 
definitions of the body axis velocities in equations (1-12) to (1-14) allows the earth-relative velocities to be 
expressed in terms of V, a, and ft: 

h = V(cos a cos ft sin 9 - sin (3 sin <f> cos 9 - sin a cos (3 cos 4> cos 9) (1-72) 

x = V [ cos a cos /3 cos 9 cos ip + sin ft (sin <j> sin 0 cos ip - cos <j> sin ip) 

+ sin a cos/?(cos<£sin 9 cos ip + sin <p sin ip)] (l' 73 ) 


y = l/[ cos a cos ft cos 6 sin ip + sin ft ( cos cp cos V* T s * n V* s * n ^ VO 

+ sin a cos ft(cos 4> sin 9 sin ip - sin <p cos i’)] 


(1-74) 


1.3 Nonlinear Observation Equations 


No standard set of observation variables exists for the aircraft analysis and control design problem. However, 
for any guidance and control problem, the main observation variables generally will be a subset of the state 
variables. Other common observation variables are the vehicle body axis translational accelerations and 
air data parameters. Thus, the dimension of g[x(f), x(t), 11 (f)] is not fixed and varies from application to 
application. The set of observation variables described in this section was selected to address a wide range 
of problems. The basic composition of the observation vector y as used in this report is given by 

y = [x T i T u T y' T ] T (1-75) 


where x and x are the state vector and time derivative of the state vector described previously, u is the 
control vector, and y' is defined by 


y' = [yi T y? yf y? yf ya T yf yST 


(1-76) 


where 


iT 

y'l “ [&x,k a y, k a z , k a x a y a z a n a x,i a y,i a z t i a n,i n \ 

y' 2 = [a M Re Re' q q c q c /p . a p a Pt T T t ] T 
ya = [7 fp a h] T 
yi = [E s P S ] T 

y's - [L D N A] t 
y' G = [u v w ii v u?] T 

y~ = [a,i ft,i h ti h,i] T 

ys = [T 1 h <7s r s ] T 


(1-77) 

(1-78) 

(1-79) 

(1-80) 

(1-81) 

(1-82) 

(1-83) 

(1-84) 
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with the elements of y* being terms related to the vehicle body axis acceleration, the elements of y 2 being 
air data terms, the elements of being flightpath-related terms, the elements of y\ being terms related to 
vehicle energy, y£ being a vehicle force vector, the elements of y f 6 being body axis translational rates and the 
time derivatives of those terms, y f 7 being a vector of variables representing measurements from instruments 
not located at the vehicle center of gravity, and the elements of y§ being a collection of miscellaneous terms. 
Obviously, this grouping of terms is somewhat arbitrary and is done primarily to ease the definition of these 
terms in the following sections of this report. This grouping of observation variables parallels that used by 
Duke and others (1987). 

The vector function g relating the state vector, the time derivative of the state vector, and the control 
vector to the observation vector is an ^-dimensional function composed of four subfunctions: 

g[x(0,x(0,u(f)] = [x T x T u T g /T ] (1-85) 

where x, x, and u are identity functions on the state vector, time derivative of the state vector, and control 
vector, respectively, and g' is composed of vector subfunctions defining the y' vector. 

The state vector, time derivative of state vector, and control vector components of the observation 
vector are not discussed in detail in this section of the report. The equations for the elements of the time 
derivative of the state vector were developed in section 1.1. The observation equations for the state and 
control variables are simply identities. The equations for the remaining observation variables are obtained 
from a variety of sources. In addition to the previously cited sources, Clancy (1975), Dommasch and 
others (1967), Gainer and Iloffman (1972), and Gracey (1980) provide the background and derivation of 
the observation equations used in this report. 


1.3.1 Accelerations. — The vehicle body axis accelerations and accelerometer outputs constitute the 

set of observation variables that, after the state variables themselves, are most important in the aircraft 
control analysis and design problem. These accelerations and accelerometer outputs are measured in units 
of g and are derived directly from the body axis forces defined in section 1.2.2. The body axis acceleration 
vector a can be expressed as 

* = T, v= i v + axr t 1 ' 86 ) 

It is important to note here that the u, u, and w body axis velocity rates, derived in appendix B and 
defined by equation (B-l), are not the body axis accelerations. The body axis accelerations contain not 
only the body axis velocity rates but also the rotational velocity and translational velocity cross-product 
terms. Thus, expanding equation (1-86) yields 





u + qw — rv 

a = 

Gy, k 

= 

v + ru — pw 


_ a>z , k 


w + pv — qu 


(1-87) 


where a y ( k, and a z y are the kinematic accelerations in the vehicle body a;, y, and 2 axes, respectively. 
Using 


( 1 - 88 ) 


ii 


(l/m)(A"T + X a + X g ) + rv - qw 

V 

= 

+ Y a + Yg) + pw - ru 

w 


_ (l/m)(Z T + Z a + Z % ) + qu - pv 


(stated as eq. (B-l) in app. B), equation (1-87) can be rewritten as 


a x,\n 


' (l/m)(X T + X a + X g ) ' 

a y,k 


(i/m)(r T + y a + y g ) 

k 


(l/m)(ZT -j- Z a + Z s ) 


(1-89) 
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where X a , Y a , and Z a are total aerodynamic forces and X g , Y g , and Z g are total gravitational forces along 
the x, y, and z body axes, respectively. This can be expanded in terms of the gravitational and aerodynamic 
forces to give (in units of g) 

Xj - D cos a + L sin a — gm sin 6 

Yi + Y + gmsm<f>cos9 (1-90) 

Zj - D sin a - L cos a + gm cos <f> cos 6 

where g 0 is the acceleration due to gravity at sea level. 

The outputs of body axis accelerometers at the vehicle center of gravity are simply the body axis 
accelerations due to the thrust and aerodynamic forces. The accelerometer output equations can be written 
directly from equation (1-90) as 

[ Xt - D cos a + L sin a 

— Yr + Y (1-91) 

9o m Z T — D sin a - L cos a 

where a x , a y , and a z are the outputs of accelerometers at the vehicle center of gravity and aligned with the 
vehicle body x, y, and z axes, respectively. Because the normal acceleration a n is defined by 

a n = -a z (1-92) 

an expression for this variable can be extracted from equation (1-91): 

a n = ( — Zt + Ds'ma + Lcosa)/g 0 m (1-93) 

The equations defining the output of accelerometers aligned with the vehicle body axes but displaced 
from the vehicle center of gravity are derived by Gainer and Hoffman (19<2) using the definition of ineitial 
acceleration given in equation (1-86) 

a = + a X V 

0 1 

and the definition of inertial velocity 

V = i-r + fi X r (1-94) 

ct 

The results from Gainer and Hoffman (1972) are reproduced here without rederivation: 

a x ,i 1 a* - [(q 2 + r2 ) x x ~ (. pq ~ r)Vx ~ (pr + q)z x ]/go 

dy yl = a y + [(pq + r)x y - ( p 2 + r 2 )y y - (qr - p)z y ]/g 0 (1-95) 

a z ,i J a z + [(pr - q)x z + (qr + p)y z - ( q 2 + p 2 )z z ]/g 0 

where a Xi ,-, a y< i, and a 2?t are outputs at accelerometers aligned with the x, y, and 2 body axes but not located 
at the vehicle center of gravity; the subscripts x, y, and 2 refer to the x, y , and 2 body axes, respectively; 
and the symbols x, y, and 2 refer to the x, y, and 2 body axis locations of the sensors relative to the vehicle 
center of gravity. Because the normal acceleration is the negative of the 2 body axis accelerometer, the 
output of a normal accelerometer not at the vehicle center of gravity but aligned with the 2 body axis, 
is given by 

a n<i = a n - [(pr - q)x z + (qr + p)y z - ( q 2 + p 2 )z z ]/go (1-96) 

The final quantity included in the general category of accelerations is load factor n. This quantity is 
defined without inclusion of the 2 body axis force component as 

(1-97) 




_L_ 

mg 
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1.3.2 Air data parameters. — The air data parameters having the greatest application to aircraft 
dynamics and control problems are the sensed parameters and the reference and scaling parameters. Chosen 
for inclusion as the sensed parameters are impact pressure q c , static or free-stream pressure p a , total pressure 
p t , ambient or free-stream temperature T, and total temperature T t . The selected reference and scaling 
parameters are Mach number M, dynamic pressure q , speed of sound a, Reynolds number Re, Reynolds 
number per unit length Re', and the Mach meter calibration ratio q c /p& • The derivation of these quantities 
is treated extensively by Gracey (1980). 

The nonlinear equations defining these quantities are 


a 


M 

Re 


Re' 


q 


<lc 


9c_ 

Pa 


l PoTo J 
V 

a 

pVl 


pV 

P 



f [(1.0 + 0.2M 2 ) 3,5 - l.Ojpa (M < 1.0) 

( {1.2M 2 [5.76M 2 /(5.6M 2 - 0.8)] 2 * 5 - 1.0}p a (M > 1.0) 

f (1.0 + 0.2M 2 ) 3 5 - 1.0 (M < 1.0) 

| 1.2M 2 [5.76M 2 /(5.6M 2 - 0.8)] 2 - 5 - 1.0 (M > 1.0) 


T t = T(1.0 + 0.2M 2 ) 


(1-98) 

(1-99) 

(1-100) 

( 1 - 101 ) 

( 1 - 102 ) 

(1-103) 

(1-104) 

(1-105) 


where p is the density of the air, p is the coefficient of viscosity, and the subscript 0 refers to sea level, 
standard day conditions. Free-stream pressure, free-stream temperature, and the coefficient of viscosity are 
properties of the atmosphere and are assumed to be functions of altitude alone. 


1.3.3 Flight path- related parameters. — Included in the observation variables are what might best 
be termed flightpath-related parameters for lack of better nomenclature. These terms include flightpath 
angle 7, flightpath acceleration fpa, and vertical acceleration h . The variables are defined by the following 
equations: 


7 — sin 



(1-106) 


fpa = 


£ 

9 


(1-107) 


h — a x sin 9 — a y ^ sin (j) cos 9 — a z y cos 0 cos 9 


(1-108) 
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1.3.4 Energy-related parameters. — Two energy- related parameters are included with the observa- 
tion variables considered in this report: specific energy E Si and specific power P s , defined as 


y 2 

E s — h + — — 


P* = 


dE s 

dt 


— h + 


VV 


(1-109) 

(1-110) 


1.3.5 Force parameters. — The set of observation variables being considered also includes four force 
parameters. These quantities are total aerodynamic lift £, total aerodynamic drag D y total aerodynamic 
normal force JV, and total aerodynamic axial force A , defined as 


L = qSC L 

(1-111) 

D = qSC D 

(1-112) 

N = L cos a + D sin a 

(1-113) 

A = —L sin a + D cos a 

(1-114) 


where S is the surface area of the wing, Cl coefficient of lift, and Co coefficient of drag. 


1.3.6 Body axis rates and accelerations. — Because they are of interest in the control analysis and 
design problem, six body axis rates and accelerations are included as observation variables. These include 
the x body axis rate u> the y body axis rate u, and the z body axis rate w. Also included are the time 
derivatives of these quantities, it , u, and th, respectively. 

The definitions of the body axis rates are given in equations (1-9) to (1-11) as 


u = V cos a cos j3 
v — V sin (3 
w = V sin a cos / 3 

The time derivatives of these terms can be defined using equation (B-l) and equations (B-8), (B-9), (B-10), 
and (1-56) as 


u — 


Xt — gtn sin 9 — D cos a + L sin a 


m 


+ rV sin / 3 — qV sin a cos j3 


Yt + 9 m 4 > cos 8 + Y 


m 


+ pV sin a cos j3 — rV cos a cos / 3 


w — 


Zt + grrt cos cos 6 — D sin a — L cos a 


m 


+ qV cos a cos j3 — pV sin (3 


(1-115) 

( 1 - 116 ) 

(1-117) 


1.3.7 Instruments displaced from the vehicle center of gravity. — The need to include measure- 
ments from instruments displaced from the vehicle center of gravity arises from the fact that not all aircraft 
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instrumentation is located at the vehicle center of gravity. The most important of these quantities are un- 
doubtedly the accelerometer outputs treated in section 1.3.1. In this section four additional parameters are 
presented: angle of attack (a |t *), angle of sideslip (/3 ,), altitude (/*,;)> and altitude rate ( h it ) measurements 
from instruments displaced from center of gravity by some x y y , and z body axis distances. The subscripts 
a, /?, /i', and h refer to the displacements of the angle- of- at tack, an gle-of- sideslip, altitude, and altitude rate 
instruments from the vehicle center of gravity. The equations used to compute these quantities are 


o t i = a + 


(JX a Pilot 

V 


(1-118) 


I,. = P+ TXP y PZP (MIS) 

ft,' = h + Xh sin 0 — sin <p cos 6 — Zh cos <f> cos 6 (1-120) 

h } i == h + cos 9 + V'h sin <p sin 9 + z- h cos <f> sin 9) — <j>{y‘ h cos <p cos 6 — z- h sin <j> cos 9) (1-121) 


1.3.8 Miscellaneous observation parameters. — The final set of observation parameters considered 
in this report is a miscellaneous collection of parameters of interest in analysis and design problems. These 
parameters are total angular momentum T, stability axis roll ratep s , stability axis pitch rate q s , and stability 


axis yaw rate r s . The equations used to define these quantities are 

T = \(^V 2 - Zlxypq - 2 I xz pr -f I y q 2 - 2 I yz qr + I z r 2 ) (1-122) 

p s = p cos a + r sin a (1-123) 

& = q (i-i24) 

r s = — psina + rcosa (1-125) 


2 LINEAR SYSTEM EQUATIONS 

Tlie standard state equation for a linear differential system has the form 

x(!) = /x(f) + 5'u(f) (2-1) 

where, for a time-invariant system, A' is a constant n X n matrix and B' is a constant n X k matrix. The 
standard output equation has the form 

y(<) = H'x(t) + F'u(0 (2-2) 

where //' is a constant i x n matrix and F' is a constant £ X k matrix. The generalized linear system 
equations used with an extended formulation compatible with the generalized nonlinear equations (1-3) and 
(1-4) can be characterized by 

Cx(f) = Ax(f) + Bu(t) (2-3) 

y (0 = Hx(i ) + Gx(t) + Fu(t) (2-4) 

where C and A are constant n x n matrices, B is a constant n x k matrix, H and G are constant £ X n 
matrices, and F is a constant £xk matrix. The nonlinear system equations developed in section 1 (eqs. (1-1) 
to (1-4)) can be linearized about a trajectory, and a linear model can be formulated that is similar to either 
the standard or the generalized linear system equations. 
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2.1 Linearization of the State Equation 


If uo (t) is given input to a system described by the state differential equation (1-3), and if xo(t) is a known 
solution of the state differential equation, then approximations to the neighboring solutions can be found 
for small deviations in the initial state and in the input by using a linear state differential equation. The 
nonlinear state differential equation (1-3) can be linearized about a general trajectory, as by Kwakernaak 
and Sivan (1972) and Dieudonne (1978), so that x 0 (i) satisfies 

Tx 0 (t) = f[x 0 (t),xo(t),u 0 (t)] 

Assuming that the system is operated at close to nominal conditions with u(t), x(i), and x(t) deviating 
only slightly from u 0 (<), x 0 (f), and x 0 (<), the following expressions can be written: 

u (t) = u 0 (t) + Su(t ) (2-5) 

x(f) = x 0 (t) + Sx(t) (2-6) 

x(t) = x 0 (0 + 6x(t) (2-7) 


where 6x(t), and 6x(t) are small perturbations to the control, state, and time derivative of the state 

vectors, respectively. 

Substituting equations (2-5) to (2-7) into the nonlinear state differential equation (1-3), expanding in a 
Taylor series about x 0 (t), Xo(t), u 0 (t ), and assuming T constant with respect to x(/) yields 

T[x 0 (t) + Sx(t )] = f[x 0 (t),x 0 {t), u (t)} +—6x+—6x+—6u + h (t) (2-8) 

where df/dx, di/dx, and c?f/c?u are defined in equations (2-9) to (2-11) and h(i) represents the sum of 
the higher order terms in the Taylor series, assumed to be small with respect to the perturbations. The 
matrices used in the Taylor series expansion are defined by the following relationships: 


5f _ df_ 
dx ~ dx 

df_ _ df 
dx = dx 

df_ = df_ 

fa ~ du 

the («, J )th elements of which are defined as 




(2-9) 

(Xo,Xo,Uo) 

(2-10) 

(Xo,Xo,Uo) 

j 

(2-11) 

(Xo.Xo.Uo) 


. dt l 

dXj 

(2-12) 

■ «* ^ 
tt 

(2-13) 

II 

OD cv, 

Jr -sr* 

(2-14) 


respectively, where f; is the zth simultaneous equation of the nonlinear state differential function in equa- 
tion (1-3), Xj the jth element of the state vector, Xj the jth element of the time derivative of the state 
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vector, Uj the jftli element of the control vector, and all derivatives are evaluated at the nominal condition 

(x 0 (/), x o (0, u 0 (0)- 

Subtracting equation (1-3) from (2-8), rearranging terms and neglecting the higher order terms yields a 
linearized state equation, 

T - Sx(t) = ^ Sx(t) + «u(0 (2-15) 

where the arguments of the matrix functions have been dropped to simplify the notation and where it is 
understood that the matrices are to be evaluated along the nominal trajectory. 

Letting 


df 

c = T -ai 

(2-16) 

df 


A ~ fa 

(2-17) * 

df 

. 

I 1 3 

» 

II 

(2-18) : 

equation (2-15) can be written as 


C 6x(t) = A Sx(t) + B 6u(t) 

(2-19) ‘ 

which is precisely the formulation of the generalized state equation desired. 

- 

Prcmultiplying both sides of equation (2-19) by C 1 results in the standard form of the linearized state 

differential equation, 

Sx(t) = C 1 A6x(t) + C 1 Z?^u(<) 

(2-20) 

Letting 

“ 

A' = C~ l A 

(2-21) 

B' = C~ 1 B 

(2-22) 

equation (2-20) can be written in the more usual notation 

<5x(/) = A* 6x(t) + B f Su (t) 

(2-23) 


2.2 Linearization of the Observation Equation 


The technique used in section 2.1 to linearize the state equations can be applied to the nonlinear observation 
equation (1-4), 

y (0 = g[x(<),x(<),u(f)] 

Performing a Taylor series expansion about the nominal trajectory (xo(£)> xq(2)> u o( 0) yields 

yo(0 + Sy(t) = g[xo(t),x 0 (t),u 0 (i)] + Sx + Sx + ^ <$u -f li(£) (2-24) 

where 


dg _ dg 
dx ~ dx 


(Xo.Xo.Uq) 


(2-25) 
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Og _ <9g 
dk ~ dk 


dg _ dg 
du ~ du 


(X 0 ,x 0) u 0 ) 


(X 0 l x 0 ,u 0 ) 


the (i, j)th elements of which are defined by 


<3g\ _ dgi 

dxJij d*j 


dkJij dkj 


(2-26) 


(2-27) 


(2-28) 


(2-29) 


(P) = |S- (2-30) 

\duJij duj 

respectively, where g; is the ith simultaneous equation of the nonlinear observation equation (1-4). Again, 
all derivatives are evaluated at the nominal condition (xo(t), ko(t), uo(t)). 

Subtracting equation (1-4) from equation (2-24), rearranging terms, and neglecting higher order terms 
results in a linear observation equation, 


where the arguments of the matrix functions have been dropped to simplify notation. Letting 

rr 


Q — 

G ~ dk 


equation (2-31) can be rewritten as 


(2-31) 


(2-32) 


(2-33) 


(2-34) 


6y(t) = II Sx(t) + G Sx(t) + F Su(t) (2-35) 

which is the generalized linear observation equation desired. 

The standard form of the observation equation can be derived by substituting for 6k from equation (2-23) 
into equation (2-33). This substitution results in 

6v(t) = II Sx(t ) + G\A' 8x(t) + D' 8u(t)] + F 6u(t) (2-36) 


which can be written as 


By letting 


8y(t) = II 8x(t) + G[A' 8x(t) + B' 8u(t)] + F 8u(t) 
8y(t) = [II + GA'\ 6x(t) + [F + GB'] 8u{t) 


II' = 11 + GA! 
F' = F + GB' 


(2-38) 

(2-39) 


equation (2-37) becomes 
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Sy(t) = II' 8x(t ) + F' 8u(t) 


(2-40) 


2.3 Definition of Matrices in Linearized System Equations 


The results of sections 2.1 and 2.2 can be used to define the matrices in the linearized system equations 
in terms of partial derivatives of the nonlinear state and observation functions taken with respect to the 
state, time derivative of state, and control vectors. All derivatives are understood to be evaluated along the 
nominal trajectory. 

Using the nonlinear state equation (1-3), 


Tx(t) = f[x(i),x(*),u(<)] 


the terms in the generalized form of the linearized state equation (2-19), 


can be defined as 


C 8x(t) = A 8x(t) + B <5u(<) 


C = 


T - 


di 

dx 


A = 


di 

dx 


B = 


di 

0 u 


(2-41) 

(2-42) 

(2-43) 


The terms in the standard form of the linearized state equation (2-20), 

8x(t) = A' 6x(t) + B' 6u(t ) 

can be defined as 



r er. 

~ 1 di 


A f = 

T — .. r 
dx J 

dx 

(2-44) 


r di] 

- 1 di 


B f = 

dx\ 

5u 

(2-45) 


In a similar manner, the nonlinear observation equation (1-4), 


y (0 = g[x(0,x(0,u(*)] 


can be used to define the terms of the generalized linearized observation equation (2-35), 

<5y(<) = II 8x(t) + G 8x(t) + F 6u(t) 
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(2-46) 

ax 

n = 5S. 

(2-47) 

Ox 

Ou 

(2-48) 


The terms in the standard form of the linearized observation equation (2-40), 

Sy (0 = II' *x(0 + F' 6u{t) 


can be defined as 


//'=& + ® S 

Ox Ox 


Ofl -1 Of 
Ox] 9x 


T' = 


% 0g ‘ 
Ou Ox . 


T 


on - 1 of 

ax. au 


(2-49) 

(2-50) 


2.4 Elements of the Linearized System Matrices 


The elements of the linearized system matrices derived in sections 2.1 and 2.2 arc determined by applying the 
linearization method employed with the vector equations in those sections to the individual scalar equations 
constituting the vector equations that define the time derivatives of the state and observation variables. 
Thus, for a matrix, such as the state matrix A defined by equation (2-42), 


A = 


Of 

Ox 


the element occupying the ith row and jth column of A, (A);j, can be represented as 


1 j,J “ 8xj 


(2-51) 


where f, is the scalar function defining the time derivative of the ith state and Xj is the jth state. The 
individual terms used in the A, B , C, II, G, and F matrices are defined in appendix D based on the 
generalized derivatives derived in appendix C. 

Using the state vector x defined in (1-7) as 

x — [p q r V a ft <j> 6 h x i/] T 


the elements of the A matrix can be expressed as 

~d(p')/dp d(p')/dq ■ ■ ■ d(p')/dy 
d(q')/dp d(q')/dq • • ■ d(q')/dy 
A = : : : 

d(x)/dp d(x)/dq ••• d(x)/dy 
_d(y)/dp d(y)/dq ••• d(y)/dy_ 


(2-52) 
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Substituting for these partial derivatives using the terms in appendix D gives 

(l/Ix)[(qSb 2 /2V 0 )C tp + dL T /dp ( l/I x )[(qSbc/2V 0 )C (q + dL T /d q + I xzPo 

— IxyTo + Ixzqo] +2I yz q 0 + r o(Jy ~ Iz)] 

( l/Iy)[(qSbc/2V 0 )C mp + dM T /dp (l/Iy)(qSc 2 /2Vo)C mq + dM T /dq 

- 2 IxzPO - lyzQO + r Q {I z - I x )] +IxyTo ~ IyzPO ] 


The elements of the B, C, II, G , and F matrices can be determined in a similar fashion, although some 
caie must be taken in determining the elements of the matrices for the observation equation and the C : 
matrix. 

To determine the elements of the matrices for the observation equation, one must consider the definition 
of the nonlinear vector function g defining the observation variables (eq. (1-85)), 

g[x(f), x(f), u(f)] = [x T x T u T g' T ] 1 




and the definitions of the matrices for the generalized linear observation equations (2-46) to (2-48) , 

rr d S 




These matrices may be expressed using a partitioning based on the vector subfunctions of g as 



(2-54) 


(2-55) 


(2-56) 


28 



which become 


H = 


G = 


1-12x12 
0l2xl2 
0fcxl2 

f* . 

0l2xl2 
Il2xl2 
0jtxl2 

. ax J 
0l2x^c 

0 1 2 X /c 

Uxfc 

L’i"j 

upon evaluating the partial derivatives of the identity functions x, x, and u. 
The C matrix may be viewed as a partitioned matrix as 


F = 


(2-57) 


(2-58) 


(2-59) 


Cn I 

l_ 

C 12 

i 

1 

i 

1 

1 

<0 

X 

<£> 

0 

1 

03X3 | 

C 22 

! 

i 

1 

_ L _ 

0gx6 

1 

1 

1 

1 

1 

1 

lcx6 

L 


1 

i 

- 


where, from equation (1-48), 


C11 


and 


J = 


1.0 ~Ixy! Ix ^xzl^x 
~ Ixz/Iy 1-0 —• Ixz/Iy 

-Ixzth -lyzlh 1.0 


(2-60) 


(2-61) 


Cn 


-d{jf)/dV —d(p')/da -d(p')/d0 
-d(q')/dV —d(q')/da -d(q')/dp 
-d(r>)/d V -d(r')/da -d(r')/d/3 


■0 ~(qSbc/2V 0 I x )Ce a ~(qSb 2 /2V 0 T x )C^- 
0 -(qSc 2 /2V 0 I y )C m . -(qSbc/2V 0 Iy)C m . 

0 -(qSbc/2V 0 T z )C n . -(qSb 2 /2V 0 T z )C n - 


(2-62) 
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C 22 = 


1.0 - 3(V)/3V -d(V)/da -d(V)/30 
—d (a)/3V 1.0 -d(a)/da -d(a)/dj3 
~30)/OV - 30)1 da 1.0 - 3(1 3)/30 m 

1.0 (qSc/2Vom)(cos (3 q C Da - sin/3 0 Cy a ) (qSb/2V 0 m)(cos j3 0 CpA 


0 1.0 + (q S c/ 2V^ m cos (3 0 )C l ^ 


(qSb/2 Vq m cos /? 0 ) Cl , 


0 (qSc/2Vgm)(sm/3 0 Cjd a + cos/? 0 C Yd ) 1.0 - (qSb/2V$m)(sm(5 Q Cn + cos/3 0 Cy.) J 


(2-63) 


The inverse of the C matrix, C x , can be expressed as a partitioned matrix in terms of the matrix subpar- 
titions of the C matrix as 



— 



-C£C l2 C£ 

1 

1 

i 

06X6 


1 _ 

03X3 j 


1 

1 

i 

_U_ 




0 6x6 

I 

l 

i 

1 

l 

lex6 



_ 


i 


_ 


(2-64) 


The elements of the A\ B\ //', and F? matrices can be determined using the C~ l matrix defined in 
equation (2-64), the A, B , 77, £7, and F matrices, and the definitions for A\ B f , 77', and F ' given in 
equations (2-21), (2-22), (2-38), and (2-39). 

3 CONCLUDING REMARKS 

This report derives and defines a set of linearized system matrices for a rigid aircraft of constant mass, flying 
in a stationary atmosphere over a flat, nonrotating earth. Both generalized and standard linear system 
equations are derived from nonlinear six-degree-of-freedom equations of motion and a large collection of 
nonlinear observation (measurement) equations. 

This derivation of a linear model is general and makes no assumptions on either the reference (nominal) 
trajectory about which the model is linearized or the symmetry of the vehicle mass and aerodynamic 
properties. 


Arnes Research Center 

Dry den Flight Research Facility 

National Aeronautics and Space Administration 

Edwards, California , January 8, 19S7 
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APPENDIX A — AERODYNAMIC FORCES AND MOMENTS 


The aerodynamic forces and moments acting on an aircraft are the result of multiple factors whose signif- 
icance varies with flight condition as well as from vehicle to vehicle. In general, these forces and moments 
are nonlinear functions primarily of Mach number, angle of attach, angle of sideslip, altitude, rotational 
rates, and control-surface deflections. For the purposes of this report, the aerodynamic forces and moments 
are assumed to be functions having the following form: 

F = $(a,P,V,h,p,q,r,a,/3,6i,...,6n) (A-l) 

where F is an arbitrary force or moment, $ is an arbitrary function, and the Si are the n control surface 
deflections. These forces and moments are related to the nondimensional force and moment coefficients by 
the equations for the forces, 

D = qSC D 
Y — qSCy 
L = qSCh 

and the moments, 

L = qSbCt 
M = qScCm 
N = qSbC n 

where b is reference span and c is reference aerodynamic chord. 

While the nondimensional aerodynamic force and moment coefficients are themselves nonlinear func- 
tions of the vehicle states, time derivatives of the vehicle states, and the control surface deflections, these 
coefficients are commonly expressed in linear form in terms of partial derivatives of these coefficients with 
respect to the functional variables. These linear equations for the aerodynamic force and moment coeffi- 
cients are derived in the same way as the linearized system equations (section 2); therefore, this derivation 
will not be repeated here. These linear equations are 

Cl = Cl 0 + Cl q « + Cl^/3 + C^Ji + Cl v V 

+ £ Cl s . A + C l p p + C L J + C Lr f + C La d + c H b (A-S) 

i=l 

Cd = C Do + Co a o + Cd^/ 3 + Co h h + Cd v F 

+ £CD i .<5, + CD p p + CD q q + C Dr r + C D(i d + CD^/? (A-9) 

*=i 

Cy = Cy 0 -(- Cy a a + Cy 0 (3 + Cyji + Cy v V 

+ £ Cy^. Si + Cy p p + Cy q q + Cy r f + C Yq « -(- Cy (A-10) 

1 = 1 

Cl = C (o + C fa a + Ct'P + c (h h + C (v V 

+ £ C tt . S l + c (p p + C tq g + C tr f + C e J + Ctjb ( A- 1 1 ) 

1 = 1 


(A-2) 

(A-3) 

(A-l) 


(A-5) 
(A-6) 
(A- 7) 
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(A- 12) 


Cm = C mo + Cm a a + CmpP + C mh h + C my V 

n - 

+ ^ ^ Cm 6 . + Cm p p + Cm q 4 "I" C rrir V -(- Cm^Sx + 

i'=l 

C n = C n Q + C na a + C np /3 + C Uh h + C ny V 

+ ^ C n6 . Si + C np P + C nq q + C Ur T + Cn^Ot + C n ^ (A-13) 

1 = 1 

where Q 0 is the value of the coefficient along the nominal trajectory and the notation C( x is defined as 

(A- 14) 


r _ dc i 


dx 

with C i being an arbitrary force or moment coefficient and x being an arbitrary state, time derivative of state, 
or control-related parameter that for the usual derivatives is nondimensionah However, the derivatives with 
respect to altitude and velocity are not taken with respect to a nondimensional quantity. The definitions 
of these nondimensional stability and control derivatives are given in terms of the coefficient Q. The 
nondimensional stability derivatives are defined as 



dCz 

da 

(A 15) 

K 

Cit = 

OC'z 

d/3 

(A 16) 

* 


OCz 

(A 17) 

s 

dibpj 2V 0 ) 


Q, = 

0C i 

(A- 18) 

- 

d(cqj 2V 0 ) 


CfrS 

dCi 

(A- 19) 


d(br/2V 0 ) 


Qd = 

dCi 

(A-20) 

- 

0(ca/2V 0 ) 



c>Q 

(A-21) 

— 

d(b0/ 2V 0 ) 



The two other stability derivatives are not nondimensional and are defined as 

Qv = 

* 


The control derivatives are defined as 


Q, = 


0C(_ 

dV 

(A-22) 

OCz 

Oh 

(A-23) 

OCt 

(A-24) 
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The rotational terms in equations (A-8) to (A- 13) are nondimcnsional versions of the corresponding vari- 
able with 


P = 
<7 = 
r = 
a 


bp 

2Vb 

eg 

2 Vo 
br 

ZVo 

ca 

Wo 


(A-25) 
(A-2C) 
(A- 2 7) 
(A-28) 


P 


bP_ 

2F 0 


(A-29) 


Because the C(_ 0 terms are included, the force and moment coefficients are total force and moment coefficients. 
The state, time derivative of state, and control parameters on the right-hand side of equations (A-8) to 
(A- 13) are differentials. 
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APPENDIX B— DERIVATION OF THE WIND AXIS 
TRANSLATIONAL PARAMETERS V, a, AND @ 

The derivation of the wind axis translational acceleration parameters is based primarily on the definitions 
in equations (1-9) to (1-14), the body axis translational acceleration equations (1-56), and the expression 
of the force terms defined in equation (1-53). In the following sections, each of the wind axis transla- 
tional acceleration terms is derived separately after stating some preliminary definitions applicable to all 
calculations. 


B.l Preliminary Definitions 


Equation (1-56), 


— V = — F - ft x V 
St m 


can be expanded, using equations (1-54), (1-55), and (1-26), to 


u 


’(1 /m)(Xr + X a + Xg) + rv - qw 

V 

= 

(i/7n)(y T + y a + y g ) + pw - m 

w 


(l/m)(Z T + y a + Z s ) + qu - pv 


The body axis aerodynamic forces can be rewritten in terms of the stability axis forces lift Z, drag D , and 
sideforce Y : 


X a = —D cos a + Z sin a 
Pa - Y 

Z a = — D sin a — Z cos a 

The gravitational forces can be resolved into body axis components such that 


(B-2) 

(B-3) 

(B-4) 


X g — —mg sin 9 
1 g = mg sin <f> cos 8 
Z g = mg cos cf> cos 6 

These equations will be used in the derivations of the F, d, and /? equations. Thus, 
body axes can be defined and expanded as 

EA" = Xx - D cos a + L sin a - gm sin 8 

EY = Yt + Y + gm sin <f> cos 9 

EZ = Zx — D sin a — L cos a + gm cos (f> cos 9 


(B-5) 

(B-6) 

(B-7) 

the total forces in the 


(B-S) 
(B-9) 
(B- 10) 


B.2 Derivation of V Equation 

Beginning with the definition of V in terms of u, v, and w in equation (1-12), 

V = (u 2 + v 2 + w 2 ) 1/2 


PRECEDING PAGE BLANK NOT FILMED 
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the equation for V becomes 


(B-ll) 


l> = s ' / = l <“ 2+ ” 2 + ”’ 2 ) ,/2 

which after expanding the derivative and cancelling terms, becomes 

V = -l(tm + vv + ww) (B-12) 

Cy substituting the definitions for u, v, and w from equations (1-9) to (1-11) and cancelling terms, equa- 
tion (B-1'2) yields 

V = u cos q cos/3 + usin/? + w sin a cos/? (B-13) 


The definitions for u, v, and w in equation (B-l) are now used with equation (B-13) to give 

• cosacos/?. v __ 

V = — (X a + At + Ag) T cos cr cos p{rv — qw) 

+ " + Tx + Yg) + sin P(pw - ru ) 

sin a cos j3 


+ 


m 


-( A a + Z'l + Z g ) + sin a cos / 3(qu — pv) 


(B-14) 


Expanding (B-14) in terms of equations (B-2) through (B-7) and cancelling yields 

V = ~[ — D cos /? + Y sin /? + At cos a cos P + It sin P + At sin cr cos /? 

— m< 7 (cos a cos P sin 8 - sin /? sin <f> cos 6 - sin a cos /? cos 4> cos 0)] 

+ tv cos a cos P - qw cos a cos /? + pw sin ft — ru sin ft 

+ qu sin a cos /? - pr> sin a cos /? (C-15) 

Equation (B-15) can be simplified by recognizing that the terms involving the vehicle rotational rates are 
identically zero, which becomes obvious after substituting for u, v, and w in these terms. Thus, the final 
equation becomes 


V = — [- Dcosfi + Y sin/3 + X T cos a cos/? + y T sin /? + Z T sin a cos/3 
- mg(cos a cos P sin 6 - sin P sin <f> cos 6 - sin a cos p cos $ cos 0)\ 


(B-16) 


B.3 Derivation of a Equation 


The equation for a can be derived from the definition of a in equation (1-13), 

w 

a — tan — 
u 

Taking the derivative of a with respect to time, 


d 


d 


a — —a — — tan 
dt dt 


-l 


w 


u 


lien expanding and cancelling terms, the equation becomes 

1 


a = 


u 2 -f w 2 


■{uw — uw ) 


(B-17) 


(B-18) 
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(B- 19) 


Substituting the definitions of u and w from equations (1-9) and (1-11) into equation (B-18) gives 

w cos a — u sin a 
a V cos 0 

Using equation (B-l) to substitute for u and w and equations (B-8) to (B-10) to define the forces, 
equation (B-19) becomes, after rearranging terms, 

a = [-L + Z'l cos a - Xx sin a + mg(cos a cos <f> cos 9 + sin a sin 0)] 

Vm cos 0 

-| - (qu cos a — pv cos a - rv sin a + qw sin a) (B-20) 

V cos/3 

which after substituting for u, v, and w from equations (1-9) to (1-11) and combining terms gives 

a — i [ —L + Zx cos a - Xx sin a + mg( cos a cos 4> cos 9 + sin a sin 0)] 

Vm cos/3 

+ q — tan /3 (p cos a + r sin a) (B-21) 


B.4 Derivation of /? Equation 

The equation for 0 is derived from the definition of 0 as given in equation (1-14), 


sin 


V 


Taking the derivative of 0 with respect to time yields 


— sin 77 
dt V 


(B-22) 


which becomes, after expanding the derivative, substituting for V, and cancelling, 

0 = i[-u cos a sin 0 + v cos 0 — w sin a sin 0} (B-23) 

Using equation (B-l) to substitute for u, v, and w and equations (B-8) to (B-10) to define the forces, 

0 = — [ - cos a sin 0 (-D cos a + L sin a + Xx - mg sin 8) + cos/3 (Y + Yx + sin <f > cos 9) 
mV 

— sin a sin j3 {—D sin a — L cos a + + mg cos <f) cos 6 )] 


+ ^ -[ — cos a sin 0 (rv — qw) + cos 0 ( pw — ru) — sin a sin 0 (qu — pv)] 

Substituting into equation (B-24) for u, v, and w and rearranging terms yields the final equation 

0 - — !— [ D sin 0 + Y cos 0 - Xx cos a sin/3 + Yx cos 0 - Zj sin a sin 0 
mV 

+ mg ( cos a sin 0 sin 8 + cos 0 sin <j> cos 9 - sin a sin 0 cos 4> cos 0)] 

+ p sin a — t cos a 


(B-24) 


(B-25) 
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APPENDIX C — GENERALIZED DERIVATIVES 


The equations defining the time derivatives of the state variables (derived in sections 1.2.1 to 1.2.4) and those 
defining the observation variables (presented in sections 1.3.1 to 1.3.8) are used to determine the generalized 
partial derivatives of the quantities with respect to a dummy variable £. The purpose of these generalized 
derivatives is primarily to facilitate the derivation of the terms in the linearized equations presented in 
section 2.4; however, these equations have also proved to be useful for computer programs and were used 
to verify the results obtained using LINEAR (see Duke and others, 1987). 


C.l Generalized Derivatives of the Time Derivatives of State Variables 

Equations (1-39) to (1-41) define the rotational accelerations of the vehicle. These equations are used to 
determine the generalized derivatives of these quantities. 

diP) 


dZ 


1 ( r dL _ 

= tear 1 w + 2 


dM T dN _ dL T _ DM r 

+ h + h -7TT - + h ~ + '3 


dZ 


dZ 


dZ 


dZ 


ONy 

~w 


— [2p(I xz l2 - hyh) — q{Ixzh — lyzh ~ D Z I 3) + r(I xy Iy + Dyli — lyzh)] 
+ \p{hzl l — lyzl 2 — D Z 1 3) + 2q(I yz I\ — I X yl 3) — r{D x I\ — I X y I 2 + I X y I3)] 


(C-l) 


dZ 


= j-{ 

det I l 


— [p(I X yh + Dyh — Iyzl 3 ) + q(DIl — Ixyh + hzh) + 2r ( lyzh — Ixzl 2 )] 

T dL T DM T dN T dL^ T dM t ( r dNy 

h dz +l 4 ~dz +h ^z +l 2 ^r +h ~dr +h ^r 

— [2p(I X zh - Ixyh) ~ q(hzh - lyzh ~ D z h ) + r(I X yh + T) y J\ - I yz h )] 

+ \p(hzh ~ lyzh ~ Dzh) T 2q(I yz h — Ixyh) — 1'(D x h hyh T hyh)} 

dr) 

— [p(I X yh + D y h — lyzh) + q(D x h ~ hyh + hzh) + 2r[I yz h — hzh] (C-2) 


dp 

dZ 

dq 

dZ 

dr 

dZ 

Dp 

dZ 

dq 

dZ 


d(r) 1 f , dL . , dM . , dN _ r T 0 My _ T dN-y 


dZ 


= J_( 

det / \ 


h dz +h ~dz 


+ l6 ~dz +h ~dz 


+ h 


dZ 


+ h 


dZ 


dp 


— [2p(Ixzh — Ixyh) — qilxzh “ + A/^5 “ lyzh)] ^ 

+ [p(/i«73 — I yz h - D z Iq) + Mht h hyh) r(D x h Ixyh T /rz-fs)] 

5?' 1 

— [p{ hyh + Dyh — lyzh) + q(Dzh ~ Ixyh + hzh) + 2r( /y- h — /ajz/s)] ^ j (C--3) 

The quantities /j, / 2 , /j, J4, 7s, /6, /?*, Dy, T?*, and det 7 are defined in equations (1-32) to (1-38) and 
(1-42) to (1-44). 

Equation (1-50) defines the decoupled rotational accelerations of the vehicle (;/, q\ and r 1 ), which arc 
used to determine the generalized derivatives of the decoupled quantities: 


djh) 

dZ 


h 


DL DL t . r T . dp , T _ . . TX dq 

QC 4 - Q£ ~ i T hy _ qhz) + ( Phz + r h + % f lhz r h ) 

dr 1 


( Phy qly 4" 2.rl xz T qh) 


dZ\ 


(C-4) 
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1 dM dM^ , T T _ r _ dp , _ _ . 3r/ 

^ I ^ ( r ^7 + <}Tyz + 2p/ xz — rl z ) — + (r/^y — P^yz) ^ 

<3r " 

(P^x ~ 8^xy ^vlxz P^z) T^T 
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Equations (1-58) to (1-60) define the translational accelerations of the vehicle. These equations are used 
to determine the generalized derivatives of these quantities: 
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Equations (1-66) to (1-68) define the vehicle attitude rates. These equations are used to determine the 
generalized derivatives of these quantities: 
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Equations (1-72) to (1-74) define the earth-relative velocities of the vehicle. These equations are used 
to determine the generalized derivatives of these quantities: 
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C.2 Generalized Derivatives of the Observation Variables 


The vector equation (1-90) defining the body axis kinematic accelerations is used to determine the gener 
alized derivatives of the individual body axis accelerations: 
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Vector equation (1-91) defines the output of body axis accelerometers at the vehicle center of gravity 
and is used to determine the generalized derivatives of the individual body axis accelerometers: 
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Using equation (1-93), the generalized derivative of the output of a normal accelerometer at the vehicle 
center of gravity can be expressed as 
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The vector equation (1-95) defining the output of orthogonal accelerometers aligned with the body axes 
but displaced from the vehicle center of gravity is used to determine the generalized derivatives of these 
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quantities: 
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Equation (1-96) defines the output of a normal accelerometer aligned with the 2 - body axis but not located 
at vehicle center of gravity, a n> ,-. This equation is used to determine the generalized derivative of a n y: 
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In equations (C-'20) to (C-23), the partial derivatives of the vehicle rotational rates with respect to the 
dummy variable £ are defined by equations (C-l) to (C-3). The partial derivatives of the outputs of the 
body axis accelerometers at the vehicle center of gravity are defined by equations (C-16) to (C-19). In these 
equations, as before, the subscripts x, y, and z refer to the x, y, and 2 body axes, respectively, and the 
symbols x, y, and 2 refer to x, y, and z body axis locations of the sensors relative to the vehicle center 
of gravity. 

Using equation (1-97), the generalized derivative of the load factor can be defined as 


d(n) 1 dL 

d£ mg d£ 


Equations (1-98) to (1-105) define the air data parameters of interest for this report. These equations 
are used to determine the generalized derivatives of the air data parameters: 
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In tlie preceding equations, the generalized derivative of Mach number appears several times. This tcim 
can be expanded using equation (C-29). 

The definitions of the flightpath-related parameters are presented in equations (1-106) to (1-108). These 
definitions are used to derive the generalized partial derivatives of the flightpath-related parameters: 
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The partial derivatives of altitude rate h and velocity rate V that appear on the right-hand side of these 
equations are defined in equations (C-13) and (C-7), respectively. The partial derivatives of the body axis 
accelerations appearing in equation (C-38) are defined in equations (C-16) to (C-18). 
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Using equations (1-109) and (1-110), the generalized derivatives of the energy-related parameters are 
defined. The partial derivatives of altitude rate and velocity rate appearing in equation (C-40) are defined 
in equations (C-13) and (C-10), respectively: 
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The derivatives of the force parameters, lift (eq. (1-111)) and drag (eq. (1-112)), are defined in sec- 
tion D.l. The generalized derivatives of the normal force (eq. (1-113)) and the axial force (eq. (1-114)) are 
presented in terms of the generalized derivatives of the lift and drag forces: 
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The body axis rates are defined in equations (1-9) to (1-11). The time derivatives of these terms are 
defined in equations (1-115) to (1-117). These equations are used to derive the generalized derivatives of 
the body axis rates and accelerations: 
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The outputs of various instruments displaced from the vehicle center of gravity are defined in equa- 
tions (1-118) to (1-121). These equations define angle of attack, angle of sideslip, altitude, and altitude late 
instrument outputs. The generalized derivatives of the quantities are based on these equations: 
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The generalized derivatives of bank angle rate, pitch attitude rate, and altitude rate with respect to the 
dummy variable Z are defined in equations (C-10), (C-ll), and (C-13), respectively. 

The final set of observation variables is defined in equations (1-122) to (1-125). These equations, defining 
total angular momentum and the stability axis rotational rates, are used to determine the generalized 
derivatives of these quantities: 
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APPENDIX D— EVALUATION OF DERIVATIVES 

The generalized partial derivatives presented in equations (C-l) to (C-56) contain partial derivatives of the 
state variables, thrust forces, and total aerodynamic forces and moments with respect to the dummy variable 
£. In this appendix, these partial derivatives are defined with respect to specific state, time derivatives of 
state, and control variables. The derivatives of atmospheric parameters are also discussed. 

D.l Preliminary Evaluation 

First, the partial derivatives of the state variables with respect to the state, time derivatives of state, and 
control variables are considered. All partial derivatives of the state variables with respect to the state 
variables are either equal to zero or unity. Thus, 

dp dg dr dV_da_d0_d<f>_d_O__(hp_dh_d^_dy_ = ^ 
dp ~ dq ~ dr ~ dV ~ da ~ d/3 ~ d<j> ~ dO ~ dip dh dx dy 
and all other derivatives of state variables with respect to state variables are equal to zero. The partial 
derivatives of the state variables with respect to the time derivatives of the state variables (a and ;3 , in 
particular) are equal to zero. This is also true of the partial derivatives of the state variables with respect 
to the control variables. 

Second, the partial derivatives of the aerodynamic forces and moments with respect to the state, time 
derivatives of state, and control variables are evaluated. Using the definitions of the force and moment 
coefficients presented in appendix A, the partial derivatives can be explicitly evaluated in terms of the 
stability and control derivatives. 

D.1.1 Rolling moment derivatives. — 
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D.1.2 


Pitching moment derivatives. — 
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D.1.3 Yawing moment derivatives. — 
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D.1.4 


D.1.5 


Drag force derivatives. 
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Sideforce derivatives.- 
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D.1.6 Lift force derivatives. 
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Next, the partial derivatives of the powerplant-induced forces and moments with respect to the state, 
time derivative of state, and control variables are considered. The partial derivatives of the powerplant- 
induced forces and moments are assumed to be zero except for moments taken with respect to the body 
axis rates ( p , q, r), moments and forces taken with respect to the velocity and velocity orientation terms 
(V, a, P), and forces taken with respect to the control variables. These terms, assumed to be nonzero, are 
taken as primitives and not evaluated further. Thus, using F p to represent a powerplant-induced force (.Yj, 
Ft, and Z^) and M p to represent a powerplant-induced moment (Lj, My, and iVy), 
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ajid 

dF p dF p dFp dFp dM p dM p dM p dM p dM p dM p 
dv' da' dp' dt t ’ d P ' d q ’ dr ’ ov ' da ' and ~dp~ 

are taken as primitives and not evaluated further. 


The final set of partial derivatives to be discussed are the derivatives of atmospheric parameters with 
respect to the state, time derivative of state, and control variables. In this report, all atmospheric parameters 
are assumed to be functions of altitude only. Thus, except for 


dT dp dp, dp a 
dh' dh' dh' and ~dh' 

all derivatives of ambient temperature, density, viscosity, and ambient pressure are assumed to be equal 
to zero. The nonzero quantities listed previously are dependent on an atmospheric model. Clancy (1975), 
Dommasch and others (1967), Etkin (1972), and Gracey (19S0) present discussions of atmospheric models. 
Tn this report, the quantities will be taken as primitives and not evaluated further. 
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D.2 Evaluation of the Derivatives of the Time Derivatives 
of the State Variables 

The generalized derivatives of the time derivatives of the state variables arc defined in appendix C, equa- 
tions (C-l) to (C-15). In this section, these generalized derivatives are evaluated in terms of the stability 
and control derivatives, primative terms, and the state, time derivative of state, and control variables. In 
this section, the notation d(i t )fdx t is used to represent the more correct notation dfi/dxj that is employed 
in the discussion at the beginning of section 3. This notation is used because there is no convenient no- 
tation available to express these quantities clearly — particularly not the usual notation employed in flight 
mechanics texts such as Etkin (1972) and McRuer and others (1973). The notation that defines quantities 
such as L p = 0(p)/dp and M q = d(q)/dq is misleading in this context because the definitions of those terms 
(such as L p , M q ) are based on assumptions of symmetric mass distributions, symmetric aerodynamics, and 
straight and level flight, and additionally do not include derivatives with respect to atmospheric quantities. 
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D.2.2 Pitch acceleration derivatives. — 
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D.2.5 Decoupled pitch acceleration derivatives- 
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D.2.10 Roll attitude rate derivatives. — 
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D.2.14 North acceleration derivatives. — 
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D.2.15 East acceleration derivatives — 
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D.3 Evaluation of the Derivatives of the Observation Variables 

The generalized derivatives of the observation variables are defined in appendix C, in equations (C-16) to 
(C-56). In this section, these generalized derivatives are evaluated in terms of the stability and control 
derivatives, primative terms, and the state, time derivative of state, and control variables. 
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Lateral kinematic acceleration derivatives. — 
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x body axis accelerometer output derivatives. 
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y body axis accelerometer output derivatives. — 
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z body axis accelerometer output derivatives. 
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D.3.7 Normal accelerometer output derivatives. 
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D.3.8 Derivatives of x body axis accelerometer output not at the vehicle center 
of gravity. — 
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D.3.9 Derivatives of y body axis accelerometer output not at vehicle center 
of gravity. — 

d i a y,i) _ sSb _ J_ 

Op 2V 0 g 0 m CYp gV™* %%y) 


(D-393) ; 
(D-394) ; 
(D-395) ; 
(D-396) ; 
(D-397) t 
(D-398) ; 
(D-399) [ 
(D-400) . 
(D-401) ' 
(D-402) 
(D-403) : 
(d- 404) r 
(D-405) 
(D-406) 
(D-407) 
(D-408) 
(D-409) 
(D-410) 


(D-411) 


72 



d( a y,») 

dq 

d(a y ,i) 

dr 

d( a g,») 

dV 

djay/) 

da 

^(fly,t) 

dP 

d<p 

09 

fl(fly.t) 

dip 

d( a »,0 

Oh 

dx 

dy 

0(«y,«) 

dp 

d jOy^z) 

dq 

d{a y ,i) 

dr 

0{ a y>i) 

da 

0{ a y,i) 

00 

0{a y >i) 

d6; 


J 5 — Cy, + -JHPoZj, + ro^v) 

2Vo5 0 m q <7o 

0 ^ Sb C Yr + ^-{qoZy - 2r 0 y y ) 
2V Q g Q m go 

1 ( OYj 

( 5/oVoCV + qSCy v + -xrr 

?5C7y « + dW 


g 0 m 

1 

<7o m 

0 


= 0 


= 0 


= 0 


= 0 


qSC\ 0 + 


3>V\ 

fi ' dp ) 


— (jSVgCY^ + qSCYb 
g Q m \2 on 


9 o 


9 o 

gSc 

2V 0 g 0 m 

gSb 

2V 0 g 0 m ~ ' > 3 
1 




Cy, 


Po m 


qSCy 6> + 


0r T \ 

d<5, J 


(D-412) 

(D-413) 

(D-414) 

(D-415) 

(D-416) 

(D-417) 

(D-418) 

(D-419) 

(D-420) 

(D-421) 

(D-422) 

(D-423) 

(D-424) 

(D-425) 

(D-426) 

(D-427) 

(D-428) 


D.3.10 Derivatives of z body axis accelerometer output not at vehicle center 
of gravity. — 
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Derivatives of normal accelerometer output not at vehicle center 
of gravity. — 
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D.3.12 Load factor derivatives. — 
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D.3.13 Speed of sound derivatives* — - 
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D.3.14 Mach number derivatives. — 
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D.3.15 Reynolds number derivatives. — 
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D.3.18 Impact pressure derivatives. — 
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D.3.19 Mach meter calibration ratio derivatives. — 
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D.3.25 Specific power derivatives. — 
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D.3.30 z body axis rate derivatives. — 
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